Flow boundary conditions from nano- to micro-scales 
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The development of microfluidic devices has recently revived the interest in "old" problems as- 
sociated with transport at, or across, interfaces. As the characteristic sizes are decreased, the use 
of pressure gradients to transport fluids becomes problematic, and new, interface driven, methods 
must be considered. This has lead to new investigations of flow near interfaces, and to the concep- 
tion of interfaces engineered at various scales to reduce flow friction. In this review, we discuss the 
present theoretical understanding of flow past solid interfaces at different length scales. We also 
briefly discuss the corresponding phenomenon of heat transport, and the influence of surface slip on 
interface driven (e.g. electro-osmotic) flows. 
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I. INTRODUCTION 

The nature of boundary conditions of fluid at solid sur- 
faces has been revisited over the recent years. Beyond the 
fondamental understanding of the fluid-solid dynamics, 
the reason for such a strong interest lies in its potential 
applications in microfluidics. The driving of liquids in 
ever tiny channels raises a number of difficulties, one be- 
ing the huge increase in hydrodynamic resistance when 
the channel size decreases. Releasing the no-slip bound- 
ary condition at the surfaces, and thereby allowing for 
boundary slippage, would allow to bypass this stringent 
conditions by decreasing wall friction. There is therefore 
a big hope to take benefit of slippage for microfluidic 
applications. 

Many experiments have been performed on the sub- 
ject, with sometimes contradicting results. We refer to 
the review of Lauga, Brenner and Stone for an exhaustive 
discussion of the experimental approaches to investigate 
slippage While experimental investigations focused 
first on slippage on bare (atomically smooth) surfaces, 
more recent works have turned towards structured sur- 
faces, in particular superhydrophobic surfaces, character- 
ized by patterns at the micro- or nano-scales. 

In this review we concentrate on the theoretical un- 
derstanding and expectations for the fluid-solid bound- 
ary condition. In particular, we will discuss the relevant 
theoretical framework for mechanisms that take place at 
different scales, and lead to slippage at solid-fluid inter- 
faces. 



II. FLOW PAST IDEAL INTERFACES/ 

From a conceptual standpoint, the simplest situation 
that can be considered is that of a semi infinite solid, 
bounded by an atomically smooth surface, in contact 
with a liquid that occupies the second half space. By 
atomically smooth, we mean for example a dense plane 



of a perfect crystalline lattice. The location of this lattice 
plane defines the xOy plane. The theoretical question be- 
comes then simply to describe momentum transfer from 
the fluid to the solid (assumed here to be at rest) in terms 
of the velocity field existing in the fluid. Let us consider 
the case, where the fluid is undergoing a laminar, planar 
Couctte flow in the x direction, e.g. driven by a second 
flat boundary far away in the z direction. In a station- 
ary state, the stress component a xz is uniform, and the 
velocity profile of a Newtonian fluid reads, far away from 
the solid wall 



v x (z) = A {z + 



(1) 



The shear stress is 777, where 77 is the shear viscosity 
of the bulk fluid and 7 the shear rate. By definition, 
v s is the slip velocity. To achieve a complete hydrody- 
namic description of the system, a constitutive equation 
involving the interface has to be introduced. Such a con- 
stitutive equation is generally described as a "boundary 
condition" . Although we will also use this traditional vo- 
cabulary, we emphasize that it is somewhat misleading. 
In many cases, "boundary conditions" are seen as auxil- 
iary mathematical constraints that allow to solve partial 
differential equations. In our mind, the status of the 
constitutive equation relative to the interface is in every 
respect comparable to the one of the bulk equations, e.g. 

• it is a material property, involving the chemical na- 
ture of the two phases that create the ideal inter- 
face. 

• it should be amenable to the same type of mi- 
croscopic justifications and statistical mechanical 
treatments as bulk properties. 

• it can be extracted unambiguously from transport 
measurements. 

As for any constitutive equation, the writing of the 
interfacial constitutive equation involves a first step of 
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FIG. 1: Schematic representation of the definition of the slip 
length b 



modeling based on considerations of symmetry and lin- 
earity. The simplest possible relation, consistent with the 
general description for newtonian fluids, assumes that the 
tangential stress exerted on the solid surface is propor- 
tional to the slip velocity, i.e. 

<Jxz = KV S (2) 

Combining @ with the constitutive equation for the bulk 
Newtonian fluid, a xz ~ r)d z v x , one arrives at the so called 
Navier boundary condition 



-d z v x = oa z v x 

K 



(3) 



which can be used as a boundary condition (in the mathe- 
matical sense) that complements the Navier-Stokes equa- 
tion inside the fluid. This boundary condition is applied 
within one atomic distance of the first solid plane Q, 
where the stress ^ acts on the solid. By construc- 
tion, the " hydrodynamic" flow that results from using 
this boundary condition will be identical to the actual, 
macroscopic flow far away from the boundaries. At the 
microscopic level, a thin interfacial region exists in which 
the flow may differ from the hydrodynamic calculation. 
We insist that this region, in which the velocity field can 
be obtained only from molecular scale studies, is not our 
main concern here. Studies at the molecular scales may 
result into velocity profiles that are not solution of the 
Stokes equation, or need to artificially introduce a posi- 
tion dependent viscosity. 

The second equality in ([3]) introduces the slip length 
b = n/n. The significance of this length is illustrated in 
figure [TJ as the distance inside the solid to which the ve- 
locity profile has to be extrapolated to reach zero. The 
standard "no slip" boundary condition corresponds to 
6 = 0, i.e. k — ► oo . Conventionally a positive slip 
length is associated with a positive slip velocity, while 



a (less common) negative length would indicate an ap- 
parent change in the sign of the velocity field near the 
solid. We emphasize that b is not an interfacial property 
in the usual sense, rather it is the ratio of a truly in- 
terfacial property, the friction k, to a bulk property, the 
viscosity rj. In particular, b is, in principle, completely 
unrelated to the width of the interfacial region alluded 
to above. All our knowledge of interfaces between solids 
and simple liquids indicates that -if critical phenomena 
are excluded- the width of the interfacial region is at most 
of a few atomic diameters. The slip length, on the other 
hand, can (theoretically) be made to diverge by going to 
the limit of a perfectly flat (in the mathematical sense) 
solid wall. Such a procedure, on the other hand, would 
not significantly alter the width of the interfacial region. 

The determination of the interfacial constitutive re- 
lation, often described as "measurement of the slip 
length" has been the subject of numerous experimen- 
tal [H, S S 1, 1 and theoretical @, 1 J, M M M El 
investigations in the past 20 years . In spite of the ex- 
perimental difficulties associated with the definition of 
"ideal" surfaces in experiments, the combination of in- 
creasingly accurate measurement methods and of insights 
from simulation and theory has produced a rather com- 
plete understanding of the "ideal" situation discussed 
here. The results can be summarized as follows. 

(i) The linear relation ([3]) can be used to describe re- 
sults for simple liquids, within the range of experimen- 
tally accessible shear rates |7J. Numerically, non lin- 
ear behavior has been reported for shear rates of order 
10 8 s _1 for simple liquids [14j . 

(ii) The friction coefficient k, and hence the slip length 
b, depend strongly on the strength of the solid liquid 
interactions. Early measurements 0, 0] were performed 
on fluid/solid interfaces that correspond to a "wetting" 
situation, with a low contact angle. These measurements 
clearly reported a "no slip" boundary condition, with a 
slip length smaller than the molecular dimensions. 

(iii) When the interaction between liquid and solid be- 
comes weaker, so that the contact angle of the liquid on 
the substrates increases, larger values of b can be mea- 
sured, typically in the range 10-50 nanometers. 

(iv) Points (ii) and (iii) correspond to a general trend, 
but there is not a unique correlation between contact 
angle and friction. Other factors, such as the crystallinity 
of the solid, do influence the friction k pi [l5j. 

These findings can be rationalized within the context 
of a linear response theory 0, Q (see also which 
has been successfully compared to molecular simulations 

[T2|. A possible starting point is the Kubo like for- 
mula, which relates the interfacial transport coefficient 
for momentum (i.e. the friction k) to the integral of the 
autocorrelation function of the momentum flux (i.e. the 
force exerted by the liquid on the solid): 



dt(f x (t)f x (0)) 



(4) 
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where S is the surface area, and f x (i) is the tangential 
stress (force per unit area) exerted by the fluid on the 
solid at time t in the x direction [62|. Note that this 
formula can obviously be generalized to the case where re 
would be anisotropic, represented by a 2 x 2 matrix, as 
would be the case for a surface of low symmetry. We will 
limit ourselves to the standard case of isotropic friction. 
Kubo formula such as JU) cannot in general be evaluated 
analytically. The correlation function can be computed 
in direct numerical simulations, or approximated using 
the standard methods of liquid state theory. 

A simple approximation of the result of equation Q 
was derived in [9(, and allows one to quantify the main 
ingredients that enter the friction. The approximate ex- 
pression reads 

Here qo is the first wave- vector of the reciprocical lattice 
of the crystalline substrate; S^^o) is the structure factor 
of the fluid evaluated for this wave-vector, in the inter- 
facial region; Dn is a collective diffusion coefficient com- 
puted at the wave- vector qo; p{z) is the density profile of 
the fluid perpendicular to the interface, and Vfs{z) is the 
interaction potential between a molecule in the fluid and 
the solid wall. The physical content of the various factors 
entering ((SJ) is clear. The structure factor describes the 
response of the fluid to the fixed atomic corrugation of 
the solid wall. A large 5y(go) indicates " commensurabil- 
ity" of some sort between the liquid and the solid, which 
will increase momentum transfer. The diffusion constant 
sets, roughly, the time scale for the decay of stress-stress 
correlations. 

Note that when re is divided by r\ to obtain 1/6, the 
product D\\T] will be formed. For simple liquids, this 
ratio is typically equal to ksT divided by a molecular 
size cr, so that any explicit reference to the microscopic 
dynamics disappears from the resulting expression for the 
slip length. 



S\\{qo)cr dzp(z)V F s{z) 2 

It is therefore expected to be independent of the fluid 
viscosity for simple liquids. 

The last important factor is the integral, which in- 
volves the wall fluid interaction and the density profile. 
It is through this integral, that the wetting properties of 
the fluid with respect to the solid enter the friction re. As 
a crude approximation, let us assume that the density 
profile consists of a layer of density p c (the " contact den- 
sity") within a molecular thickness a, followed by a bulk 
density ps- We moreover assume van der Waals interac- 
tions with Hamaker constant Ah ■ Then the integral can 
be crudely approximated by (A H / pb<? 2 ) X (Pc/pb)- 

Gathering the different ingredients that enter equation 
([S]), it is easily seen that the only one that can display 
large variations is the "contact density", i.e. the average 



density of the liquid very close to the solid surface, typi- 
cally within one molecular size, and, to some extent, the 
ratio An/ksT. Both factors go into the same direction, 
namely large slip lengths are favored by a small liquid 
solid interaction. However, it is clear that, from a di- 
mensional viewpoint, the integral J °° dzp(z)VFs{ z ) 2 will 
always remain a " molecular" quantity, so that the length 
scale that emerges from equation Q is not expected to 
become very large unless very special state points are 
considered. In fact, our simulations show that the slip 
lengths that are achieved for simple Lennard- Jones mod- 
els vary between a few molecular diameters and 50 to 60 
molecular diameters, depending on the value of the inter- 
action parameter that defines the solid liquid attraction 
(Hamaker constant). This is in very reasonable agree- 
ment, at a qualitative level, with the most recent exper- 
imental results 0]. 

The same framework may be applied to melts of short 
polymers. In this case, the substrate wave-vector qo is 
much larger than the inverse radius of gyration of the 
polymer, so that the collective diffusion is measured at 
the monomer scale : Z>m is accordingly expected to be 
independent of the molecular weight N. This predicts 
that the ratio D^rj scales like rj and increases with the 
polymer length, similar the predictions of de Gennes for 
longer chains [17| . This prediction and the validity of Eq. 
(|S"|) has been exhaustively verified for short polymers (up 
to 16 monomers) using Molecular Dynamics simulations 
by Priezjev and Troian [12j. The experimental situation 
may be more complex since slippage will be strongly de- 
pendent on the polymer adsorption at the surface. Poly- 
mers also display strong nonlinear effects at moderate 
shear rates, which are not accounted for in our descrip- 
tion, and have been extensively studied experimentally 
[II (Til and theoretically 0. 

Another interesting example of the use of equation (j4|) 
is the recent study of the influence of electrostatic effects 
on the solid-liquid friction in ionic solutions [H [3 ■ A 
simple estimate of these effects based on Eq. ((4]) leads to 
an electrostatic friction contribution to the friction fac- 
tor, varying as re' cx t}Y?Ibi with E the surface charge (in 
units of the elementary charge e), and £b = e 2 / AireksT 
the Bjerrum length (with e the dielectric constant). In- 
terestingly this contribution is independent of the Debye 
length. 

While the theoretical and numerical descriptions con- 
verge unambiguously, a number of experimental reports 
have appeared over the years, that would appear to con- 
tradict either (i) or (ii) above, or both [l III El 13 • 
In some cases, these discrepancies can be attributed to 
experimental artefacts associated with the measurement 
apparatus [Tj, [25j. However, another possible source of 
misunderstandings and errors lies in the structuration ef- 
fects described below. The intrinsic constitutive bound- 
ary condition discussed in this section may be rather dif- 
ferent from what is probed in flow experiments at larger 
length scales. 
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We conclude with a brief discussion of a view of the 
solid liquid interface which is sometimes proposed, in 
which the interfacial region is described as a "film" of 
thickness Wi, formed by some interfacial fluid (or vapor) 
of viscosity rji- Standard "no slip" boundary conditions 
are applied at the film boundary. A straightforward cal- 
culation shows that, within this description, the friction 
k that enters the constitutive equation for the interface 
is given by 

i = (i-iU (7) 
and the equivalent slip length is 

Such a description does not have a clear microscopic 
funding. In general, it is impossible to identify a "phase" 
with distinct properties at the interface. However, it 
gives a useful insight into the sensitivity of interfacial 
transport parameters to the structure of the interface. 
For example, assuming that the viscosity % is that of a 
vapor layer, it is easily shown that a slip length of 20nm 
corresponds to an interface thickness of just 0.25 nm (for 
water), hardly the size of a single molecule. Clearly this 
"layer" cannot correspond to a real phase, but is merely 
a very schematic representation of the depletion close to 
hydrophobic solid walls. 



III. CLIMBING THE STRUCTURATION SCALE 
A. Effective slippage on patterned surfaces 

The situation of an atomically smooth surface, as de- 
scribed in the previous paragraph, involves basically two 
lengths scales : the microscopic scale, £, characterizing 
the liquid-solid interface, typically of the order of the 
atomic size; and the macroscopic size £, corresponding 
to the flow scale, e.g. the size of the channel in which the 
flow is conducted. 

Intermediate scales may however be present, for ex- 
ample when the surface exhibits large scale structura- 
tion, associated with roughness or chemical patterns on 
a scale L. Assuming a scale separation between the dif- 
ferent lengths, £ <C L < £, the microscopic information 
associated with the liquid-solid dynamics may be inte- 
grated out and summarized in the local slip length on the 
patterned surface. In this case, the hydrodynamic quan- 
tities, velocity profile and pressure, do obey the Stokes 
equation, complemented by the local Navier condition 
with local slip lengths {b} (of order £) at the interface. 

Once the solution of this set of equations has been 
calculated, it is possible to define an effective boundary 
condition, by integrating out the intermediate scale L. 
This introduces an effective slip length b c g, which may 



be defined in terms of the friction force F/ on the solid 
interface and an averaged slip velocity (V s ), as : 

F f = -A V { -p- (9) 

The averaged slip velocity (V s ) is computed in a given 
plane parallel to the interface and the effective slip length 
will thus depend on the specific choice for the location of 
this plane (e.g., top, bottom or mean of the roughness). 

This definition of b e g integrates out the modulation 
of the surface at the scale L. The effective slip length, 
6 e ff) thus defines the effective boundary condition for the 
velocity profile v, at scales larger than L : 

v s = b cS d z v x (10) 

The friction force on the surface can be obtained from 
the pressure P and viscous stress tensor a at the solid 
interface dT, : 

F f = [ (-PI + f) ■ dS (11) 

with I the unit tensor. From the previous definition, the 
calculation of the effective slip length therefore reduces 
to the resolution of the initial (pure) hydrodynamic prob- 
lem, involving the evaluation of the hydrodynamic veloc- 
ity field and pressure fields. These fields are obtained 
from Stokes equation, complemented by the local bound- 
ary condition on the interface, with local slip lengths {b}. 

As a first example of such a procedure we consider the 
simple case of flow past a wavy solid surface. For simplic- 
ity we consider a flow velocity perpendicular to the mod- 
ulation of the surface, so that the situation is effectively 
two dimensional. The effective slip length is obtained in 
terms of the geometrical parameters of the interface, the 
period A = 2?:/q and amplitude u (<C A) of the rough- 
ness : b e g = -£-f{qu). In the case where a perfect slip 
boundary condition applies locally on the surface, one 
obtains f(qu) oc (qu)~ 2 (26|, while for a local no-slip 
boundary condition on the surface, f{qu) cx (qu) 2 [27| . 
In both cases, the effective boundary condition applies 
at the averaged position of the surface. Though simple, 
these results provide some interesting insights into the 
nature of effective boundary condition : 

• The result for the perfectly slipping interface shows 
that roughness inevitably decrease slippage. The 
friction originates accordingly in the viscous dissi- 
pation as the fluid flows past the roughness pattern. 
This result was first pointed out by Richardson [28[ , 
who showed that the no-slip boundary condition 
measured at the macroscopic scale follows from the 
roughness of surfaces. Recent MD simulations con- 
firm this result [lEl]. 

• On the other hand the result for a wavy surface 
with a local no-slip boundary condition suggests 
an effective slip on a wavy no-slip surface, while an 
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increased dissipation is expected on a modulated 
surface. This is a direct consequence of the choice 
for the position of the boundary plane, where the 
effective boundary condition applies : for the above 
results the averaged position has been chosen. For 
a boundary plane position at the bottom of the in- 
terface, the effective slip length is shifted by — u and 
thus negative, in agreement with the increased dis- 
sipation on the modulated surface. This shows that 
the effective BC applies at some position between 
the bottom and top of the roughness. 

• The above results applies for a single surface. If the 
flow is conducted in a channel with two confining 
surfaces separated by a distance H, the above ex- 
pressions for the effective slip lengths are modified, 
and do depend on H ! (see eg [23] )• Such a depen- 
dence is exhibited when the gap size compares with 
the roughness period H < A. This points to the fact 
that effective boundary condition is not a charac- 
teristic of the liquid-solid interface solely, but may 
depend on the flow configuration when the surface 
modulation scale L = A becomes comparable (or 
larger) than the flow scale C = H . 



B. Slippage on superhydrophobic surfaces 

Another situation which has been the object of a strong 
interest over the recent years is that of superhydrophobic 
surfaces. These surfaces combine hydrophobicity at the 
molecular scale with roughness at intermediate scales, 
leading to the so-called "Lotus leaf effect" 30]. The 
roughness amplifies the natural non-wetting caracter of 
the surface, leading to very large contact angles - up to 
175° for a liquid drop on the surface. We consider here 
the situation where the liquid surface remains at the top 
of the roughness (Cassie state). The liquid surface is 
in contact with the solid only through a small fraction 
cj) s of the surface, while the remaining is free-standing. 
At the hydrodynamic level, this composite surface leads 
to a spatially dependent boundary condition with a no- 
slip boundary condition on the real solid-liquid interface 
(with fraction (f> s ), while the remaining surface is charac- 
terized by a perfect slip boundary condition (6 = oo) on 
the liquid- vapor interface. 

Although the hydrophobicity of the surface at the 
molecular scale tends (according to the discussion of sec- 
tion [n| to increase the slip length at small scale, the es- 
sential role of hydrophobicity is here a structural rather 
than dynamical one. The rough hydrophobic surface re- 
places part of the liquid solid interface with a stress free 
liquid vapour interface. 

This problem has been first tackled theoretically by 
Philip [3l| and more recently by Lauga and Stone [HJ and 
Cottin-Bizonne et al. [33j . Several numerical approaches 
have also been proposed either at the molecular scale, 
using molecular dynamics [34| . or at larger mesoscopic 
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FIG. 2: Schematic representation of a surface patterned with 
nanoposts, and equivalent representation in terms of flow 
boundary conditions. 



scales using finite element methods, Lattice Boltzmann 
or phase field models [35|, H|| [37|] . For a pattern composed 
of stripes, the expression for the effective slip length b c g 
depends on the direction of the flow. For a flow parallel 
to the stripes, b e g reads : 



b cfi = - log( - — — ) 

TT COS[|(1-0 S )J 



(12) 



For a flow in the direction perpendicular to the flow, the 
result is given by the above expression divided by a factor 
2. These expressions exhibit a weak (logarithmic) depen- 
dency on </>„. In practice, the logarithmic term only varies 
in the interval [0.4; 1.0] for <f> s between 0.1 and 5%, i.e. 
for feasible surfaces. 

At a qualitative level, the theoretical results for the 
stripes can be summarized by stating that the effective 
slip length b e s essentially saturates at the value fixed by 
the lateral scale of the roughness L, with an unfavorable 
prefactor (l/(2?r) in Eq. fl2]l). For example Eq. (fi~2"|) 
shows that a slip length larger than the period of the 
pattern, b c s = L, is obtained only when S < 10 -3 . A 
large slippage is therefore difficult to obtain within the 
stripe geometry. 

Another, more 'natural', geometry is that of a pattern 
of posts (see figure [2]) . In this case the hydrodynamic 
model consists in assuming a flat interface with a no-slip 
BC on the top of the posts, while the remaining of the 
surfaces obeys a shear-free BC. 

No analytical resolution of Stokes equation with this 
set of boundary conditions has been performed up to 
now. However a simple 'scaling' argument accounts for 
the effective slip leng th in the limit of a small fraction 
4> s of solid zones [3g]. The friction force on the surface 
reduces to the force on the solid zones Ff — A^sij^fg, 
with 7 S the averaged shear rate on the solid patch. To 
estimate j s , one may note that the velocity profile in the 
liquid is influenced by the solid zones only in a region 
of their size, a, in all directions : this behavior is due 



6 



to the Laplacian character of the Stokes equation obeyed 
by the fluid velocity. One therefore expects y s ~ U/a, 
where U is the slip velocity of the fluid on the free slip 
zones. Therefore F ~ A4> s rjU/a. Now the effective slip 
length is defined by Ff = AV/b c s, with V ~ U the slip 
velocity. One therefore deduces, in the limit <p s — > 0, 

&eff oc — . (13) 

For a pattern made of individual posts, <f) s — (a/L) 2 , 
and 6 c ff ~ L 2 /a ~ L/y/^. A numerical resolution of the 
Stokes equation with a pattern of no-slip square dots con- 
firms the validity of the heuristic argument above in the 
limit of small <p s (basically <f> s < 40%) (38|. It moreover 
gives the prefactor of the predicted relationship, which is 
very close to 1/tt, so that (38| (in the limit 4> s — >• 0): 

b cS ^-^= (14) 

7T V<Ps 

This simple prediction would deserve an analytical justi- 
fication, which has not been performed up to now. 

We note that the above argument is independent of 
the geometry of the pattern - posts or stripes - and one 
should recover the expression in Eq. (|12p for stripes. 
For solid zones made of stripes, we have <j) s = a/L, so 
that b Q g ~ L : up to slowly varying logarithmic terms, 
this is in good agreement with the Philip-Lauga-Stone 
prediction, b c g ~ L/ log(l/ (j> s )- As we emphasized above, 
in this geometry, the effective slip length is mainly fixed 
by the roughness period, L. 

From the experimental point of view, several exper- 
imental studies have been conducted to measure slip- 
page effects on superhydrophobic surfaces. Using pres- 
sure drop experiments, Rothstein and coworkers have re- 
ported slip lengths in the micron range on patterns of 
stripes [33, |4(| , with periodicity of the stripe in the tens 
of microns range (up to 150/im). Experiments are in good 
agreement with the prediction of hydrodynamic model 
of Philip [3l| . More recently Choi et al. performed slip 
measurements on the same geometry but with nanoscale 
patterns [4l|, showing a slip length of a few hundred of 
nanometers, again in agreement with theoretical hydro- 
dynamic predictions. 

On surfaces made of a pattern of pillars and using rhe- 
ological measurements, Choi and Kim reported slip 
lengths in the tens of micron range on a surface made of 
nanoposts, 6 e ff ~ 20//m for water. The expression for the 
slip length on dilute posts derived above, b c g ~ ^L/ \fcf Sl 
requires a very small solid fraction in their case to reach 
a 20^im slip length (with L ~ 1/zm) : <j) s ~ 0.03%. The 
contact angle on such surfaces with ultra-low solid frac- 
tions is expected to be very close to 180o: with this solid 
fraction a contact angle of 179° i sexpected, on the ba- 
sis of the Cassie relationship [43|], and assuming a con- 
tact angle of 9 = H0° on the solid posts. Contact an- 
gles ~ 180° are indeed reported by Choi and Kim for 
their 'nanoturf surface [42|]. Note that the slip length is 



strongly dependent on the solid fraction in this regime 
: for example a slip length of 4.5/im is predicted for a 
solid fraction of <fi s = 0.5%, corresponding to a contact 
angle of 175°. More recently, and using also rheological 
measurements, Truesdell et al. report huge slippage (in 
the millimeter range) for surfaces with a strip geometry 
p4| . while 4> s ~ 50% in their case. This result is at odds 
with theoretical expectations and would require further 
verifications. This points to the experimental difficul- 
ties to measure slippage effects inherent to rheological 
measurements, see Ref. [45] ]. Eventually, we mention 
the results of Joseph et al. who measured slip length on 
superhydrophobic carbon nanotube carpets using micro- 
particle velocimetry (/i-PIV)measurements. Slip lengths 
in the micron range are found and measured to be pro- 
portional to the underlying pattern period L (with fixed 
surface fraction of solid (j> s ~ 0.15 in their case) [46]. This 
results is in full agreement, both qualitative and quanti- 
tative, with theoretical predictions of Eq. (fl~4"|) . 

As an overall conclusion, these results and predictions 
do show that very large slip lengths may be obtained 
only at the expense of important efforts to obtain nano- 
engineered surfaces with very small solid fraction. Typi- 
cally, slip length larger than tens of microns are expected 
only with micro-patterned surfaces for which the mea- 
sured contact angle is larger than 178°. This is a very 
stringent condition. 

C. Effective boundary condition on porous surfaces 

In this context, it is interesting to point out the differ- 
ences and similarities between our approach of ideal or 
structured interfaces and the approach used in the field 
of porous media. Effective slip boundary conditions at 
the interface between a bulk liquid and a porous media 
have been introduced empirically by Beavers and Joseph 
E3, 51, El, and strongly resemble equation ©. They 
can be justified on the basis of matching a Brinkman 
description of the flow field in the porous media 

^e//V 2 v - VP = -|v (15) 

with the standard Stokes flow in the fluid. The resulting 
boundary condition is 

d z v x (0) = ^^(v x (0)-U) (16) 

where K is the permeability of the porous medium, and 
U = —KVP is the drainage velocity deep inside the 
porous wall. This condition is often empirically modi- 
fied to account for variations in the local permeability or 
viscosity near the interface [491 ] . From a conceptual view- 
point, this approach is different from the one described in 
section UTT as it does not attribute physical properties to 
the interface. Rather, the usual hydrodynamic approach 
is used everywhere, with particular approximations rel- 
evant to porous media. The approach described in the 
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present section is somewhat closer in spirit, with the dif- 
ference that it uses essentially exact results to describe 
the flow in the interfacial region. 

IV. ANALOGY WITH THERMAL BOUNDARY 
CONDITIONS 

An interesting analogy can be established between the 
flow of a fluid past a solid interface [50f , and the flow 
of heat (or energy) transport across an interface. Let 
us consider two media (1 and 2) separated by an ideal 
planar interface perpendicular to the z axis, and a steady 
heat flux J = Je z perpendicular to the interface. Each 
medium is assumed to be described by Fourier's law, with 
conductivity The temperature profile in medium i is 
linear 

Ti(z) = -z + Tm (17) 

The temperature jump at the interface, 7i(0) — T 2 (0), is 
the exact analog of the slip velocity. Again, note that 
Tj(0) can be defined unambiguously as extrapolations of 
the bulk temperature field in the interface region, without 
any reference to a local temperature at the interface. A 
constitutive equation for the interfacial region is written 
in analogy with equation ([2]) 

T((0)-T 2 (0) = T^Jni2 (18) 

Equation ([18]) defines the Kapitsa resistance Rk (also 
known as thermal boundary resistance) , as the transport 
coefficient describing the thermal effect of the interface. 
In order to complete the analogy with the slip length, 
a "Kapitsa length" Ik can be defined using a suitable 
thermal conductivity, e.g. k\. Ik — Rkki is then the 
the thickness of material 1 which is thermally equivalent 
to the interface. 

This transport coefficient has been extensively studied 
in the context of semiconductor physics, where the resis- 
tance can be interpreted in terms of phonon scattering 
at the interface [5l|. It is also well known in the field 
of superfluid liquid helium (sTj] (where Kapitsa originally 
introduced the notion), where the very large conductivity 
of Helium II makes the effect particularly important. A 
Kubo formula similar to Q can be formulated and reads 

[HUE! 

1 <? r°° 
Rk~ = I^tJ *< ? W«( Q )) (19) 

where q(t) is the heat flux (per unit area) across the in- 
terface at time t, and S is the interfacial area. 

For simple fluid/solid interfaces, the effect of the 
Kapitsa resistance has attracted much less attention. In 
the recent years a number of numerical and experimen- 
tal studies have appeared, motivated in particular by the 
interest in thermal properties of "nanofluids" (colloidal 



solution of oxide or metallic particles). The general trend 
from simulation results [5(J EH is that the Kapitsa length 
will depend on the interfacial or wetting properties in a 
way that is comparable to the slip length. A weaker 
affinity between the solid and the liquid, will result in 
higher resistances. Other ingredients, that are irrele- 
vant for the slip length, include the difference in acoustic 
impedance between the liquid ans the solid [54]. Exper- 
imental studies involve transient absorption or transient 
reflectivity experiments [55L [56| . Typical values found 
for G = l/R K are in the range 50 - 300GW.m- 2 .K _1 , 
with the smallest values obtained for hydrophobic inter- 
faces in contact with water [53l IBH ]. This corresponds 
to a Kapitsa length (using the heat conductivity of wa- 
ter) of order 12nm, indeed comparable to the slip length 
at similar interfaces. If an attempt to model such a re- 
sistance through an "equivalent vapor layer" is made, a 
thickness of the order of one molecular size is found. This 
shows again that such transport coefficients are intrinsi- 
cally interfacial properties, that cannot be assigned to 
the existence of a different "phase" at the interface. 

We are not aware of any theoretical or experimental 
studies of heat transport across structured interfaces sim- 
ilar to those studied above for the slip properties. Obvi- 
ously, strong effects, involving thermally driven dewet- 
ting, could be expected in the case of structured hy- 
drophobic surfaces. 



V. THE EXPECTED BENEFITS OF SLIPPAGE 
A. permeability and hydrodynamic dispersion 

One of the reason for the recent interest in slippage is 
that it may facilitate flows in micrometric chanel, e.g. for 
microfluidic purposes. A simple calculation shows that 
for a pressure drop flow in a slit with thickness h and 
slip length b on the surfaces, the mean velocity is in- 
creased by a factor 1 + 6b/h as compared to the no-slip 
surfaces |63| |. Slippage thus increases the permeability 
of channels and porous materials. A second interesting 
property of slip flows is the reduction of hydrodynamic 
dispersion of a probe transported in the channel using 
a pressure drop flow. This dispersion originates in the 
vanishing velocity at the surface, while it is maximum 
in the middle of the channel. For slipping surfaces, the 
velocity does not vanish at the surface, thereby reduc- 
ing the difference between the maximum and minimum 
velocities. A criterium for dispersion may be defined as 
Av = {v m ax - «mm)/»man which for a slipping channel 
equals Av = (1 + 46//i) -1 . 

For both phenomena, the efficiency of slippage is de- 
termined by the ratio b/h. This shows that micrometric 
slip lengths have to be reached for the slippage effects 
to be relevant for microfluidic purposes. As discussed 
above, this requires a specific engineering of the surfaces 
since slip lengths on bare surfaces are well below the mi- 
cron (b ~ 20nm for water on hydrophobic surfaces). For 
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FIG. 3: Schematic representation of electro-osmotic flow past 
a slippery surface. 



example, for a channel with size H ~ 10/im and a slip 
length b ~ 20nm, the increase of permeability is of order 
of one per cent ! For a slip length b ~ 2/im as measured 
e.g. in [461 ], this permeability is more than twice its value 
for no-slip surfaces. 



B. Interfacial transport phenomena 

The above discussion relies on the fact that the flow is 
driven on the scale of the channel h. However alternative 
methods do generate flows within the interfacial struc- 
ture at the solid- liquid interface. Electro-osmosis, i.e. 
flow generation by an electric field, is one such example 
which is commonly used in microfluidics. As discussed 
in references [13, HH, [59j] these interfacially driven trans- 
port phenomena are considerably amplified by surface 
slippage, even for nanometric slip length. The reason for 
this amplification is that slip length b now compares with 
the thickness of the interfacial structure A, and one ex- 
pects a amplification ratio for the transport of order b/X. 
For typical conditions, this ratio is larger than one even 
for nanometric slip lengths and slippage strongly increase 
the efficiency of interfacial transport phenomena. 

Let us precise the underlying mechanisms on the ex- 
ample of electro-osmosis. Under the application of an 
electric field E, the solvent acquires a plug flow like ve- 
locity profile, with a velocity veo proportional to E and 
given by the Smoluchowski formula : 



veo 



V 



(20) 



with e the dielectric permittivity of the solvent, n its 
viscosity; £ is the Zeta potential of the surface, which is 
traditionally assumed to match the electrostatic potential 
at the position where the velocity profile vanishes (the 
shear-plane position) [60l ] . 

Electro-osmosis originates within the Debye layer at 
the interface [6(3]. The Debye layer quantifies the width 
of the interfacial regions at a charged surface and results 



from the competition between the attraction of ions at 
the charged surface and entropic effects. Its typical size 
Ad depends on the salt concentration in the liquid, and 
is of the order of a few nanometers : Ad ~ 3nm for a con- 
centration of salt of 10~ 2 molar. The interfacial Debye 
layer is electrically charged, while the remaining bulk sol- 
vent is electrically neutral. Therefore, when applying an 
electric field, the electric driving force is located within 
the Debye layer only. The driving force per unit sur- 
face acting on the liquid thus writes f e = aE with a the 
charge of the Debye layer, which is exactly opposite to the 
surface charge a ~ cVo/Ad, for weak surface potentials 
Vq. This leads to f e ~ —ej^E. Now the electro-osmotic 
velocity results from the balance between this electric 
driving force in the Debye layer and the viscous stress at 
the surface. For a no-slip boundary condition, the lat- 
ter is f n ~ rjVEO /A d , since the velocity varies on a scale 
given by the Debye length. If slippage is exhibited at the 
surface, the velocity varies on the size Ad + b and the 
viscous stress becomes accordingly f v — t]Veo/(Xd + b), 
lower than for the no-slip BC. Gathering these results we 
get 



n 



veo 



Xj 



Ad 



(21) 



This simple argument leads to the Smoluchowski formula, 
Eq. (|20p , and provides an expression of the zeta potential 
in the form 



C = V 1 



b_ 
Ad" 



(22) 



This expression can be obtained using more rigourous 
description, see [2l|, [5?], . Molecular dynamics simula- 
tions of charge transport confirm fully this relationship 
[2lll57l]. They show, furthermore, that for a non-slipping 
surface Vq can be identified with the potential at the 
plane of shear Vq = V(z s ), with z s of the order of one 
liquid layer. As expected a strong amplification of charge 
transport is therefore demonstrated on slipping surfaces. 

This behaviour can be generalized to all interfacial 
transport phenomena. As practical examples, one may 
cite diffusio-osmosis and thermo-osmosis. These two phe- 
nomena correspond to the induction of a flow by the gra- 
dient of a solute concentration for the former and by a 
gradient of temperature for the latter [6l| . The flow ve- 
locity for these phenomena is found to be proportional to 
the applied gradient of the observable O under consider- 
ation (electric potential, concentration, temperature) : 



v s = --TL.[l + b/L}.^- 
7] ax 



(23) 



r and L defined in terms of the (equilibrium) inter- 
facial stress anisotropy, Y, eq {x,y) = a n — o~ t : T = 

J^dy-^WT (V)' and L = r- 1 J™dy. V .^(v). These 
quantities are length scales of the order of the range of 
interaction of the solute with the solid surface [53 
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The origin of the phenomena thus lies in the thin inter- 
face layer of size L where the solute interacts specifically 
with the surface. For the diffusio-osmosis phenomena, 
the driving force is an osmotic pressure gradient within 
the interface layer. For thermo-osmosis, this is surface 
energy gradient induced by the temperature dependence 
of the stress anisotropy. 

For all these phenomena, slippage again amplifies the 
transport phenomena with a ratio b/L. Slippage, even 
nanometric, is therefore relevant to amplify interfacial 
transport in microfluidic systems. 

VI. PERSPECTIVES 

In this short review, we have tried to establish a clear 
picture of our present understanding of slippage at differ- 
ent scales. At the molecular scale, we consider surfaces 
as ideal and introduce an " intrinsic" constitutive relation 
for the interface. The corresponding slip length depends 
strongly on the nature of intermolecular interactions (hy- 
drophobicity effect), but stays within "molecular" orders 
of magnitude, a few ten nanometers at most. 

A different picture is obtained by considering surfaces 
structured on an intermediate scale. To describe such 
situations, we adopt a mesoscopic view of the interface, 
in which the molecular details are hidden in a local, po- 
sition dependent slip length. Such a view can be val- 
idated through molecular simulations. Hydrophobicity 
effects are essential to reduce the contact area between 
solid and liquid, so that a significant reduction of the 
solid-liquid friction is obtained. When this reduction is 
achieved, much larger slip lengths can be observed. The 
slip length is essentially determined by the scale of the 
surface pattern and the surface area of the "dewetted" 
part, on which a perfect slip boundary condition applies 
at the local level. 

This description at two different levels is necessary to 
describe the full complexity of surfaces that can now be 
engineered at a submicrometer scale, with a combination 
of local hydrophobic properties and of a complex geo- 
metric pattern. For patterns whose typical scale is much 
larger than the intrinsic slip length, the separation of 



lengths scales allows a complete decoupling between the 
molecular and the hydrodynamic descriptions. Molecu- 
lar simulations indicate that the same type of approach 
is possible, even when the length scale of the pattern and 
the intrinsic slip length are similar. 

From this complete theoretical description, it is clear 
that the control of actual flow properties through surface 
engineering necessitates a control of the surface state at 
each different scale. When such a control is achieved, 
the flow of simple liquids appear to be relatively well 
understood, both experimentally and theoretically. More 
complex situations, such as the amplification of interface 
driven flows by molecular slippage described in section 
IV B[ have not been studied experimentally in detail yet, 
but are promising for the design of microfluidic devices. 

In section HVl we discussed the analogy between "tem- 
perature slip" and " velocity slip" . Experiments and sim- 
ulations have already shown that this analogy can be 
made quantitative, at least when ideal surfaces are con- 
cerned. The role of nano-structuration of the surfaces 
for thermal transport has not been explored up to now. 
Nonlinear effects may be expected for strong heat fluxes, 
while cross effects between heat and momentum trans- 
port are expected to be negligible, based on symmetry 
arguments. 

Finally, it could be interesting, at least from a formal 
point of view, to extend the above discussion of heat and 
momentum transport across an interface to diffusive mass 
transport. The standard description assumes a local par- 
tition equilibrium at the interface, or equivalently a con- 
tinuity of the chemical potential. Following the approach 
outlined here, it would seem more consistent to introduce 
a discontinuity in chemical potential proportional to the 
diffusive flux across the interface. The magnitude of this 
discontinuity has not, to our knowledge, be explored at 
the molecular scale. 
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